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We theoretically investigate the high-temperature thermodynamics of a strongly interacting 
trapped Fermi gas near either s-wave or p-wave Feshbach resonances, using a second order quan- 
tum virial expansion. The second virial coefficient is calculated based on the energy spectrum of 
two interacting fermions in a harmonic trap. We consider both isotropic and anisotropic harmonic 
potentials. For the two-fermion interaction, either s-wave or p-wave, we use a pseudopotential 
parametrized by a scattering length and an effective range. This turns out to be the simplest way 
of encoding the energy dependence of the low-energy scattering amplitude or phase shift. This 
treatment of the pseudopotential can be easily generalized to higher partial-wave interactions. We 
discuss how the second virial coefficient and thermodynamics are affected by the existence of these 
finite-range interaction effects. The virial expansion result for a strongly interacting s-wave Fermi 
gas has already been proved very useful. In the case of p-wave interactions, our results for the high- 
temperature equation of state are applicable to future high-precision thermodynamic measurements 
for a spin-polarized Fermi gas near a p-wave Feshbach resonance. 



I. INTRODUCTION 

Strongly interacting fermions occur in many fields of 
physics, ranging from condensed matter physics to nu- 
clear physics, astrophysics and cosmology. Recently, a 
new type of fermionic superfluid has been realized in ul- 
tracold atomic gases of "^^K and ^Li confined in a har- 
monic trap. In these systems the coupling strength be- 
tween the fermions can be tuned precisely by an exter- 
nal magnetic field from weak to infinitely strong across 
a collisional (Feshbach) resonance. In the case of s- 
wave interactions, this has already led to the observation 
of a smooth crossover from a Bardeen-Cooper-Schricffer 
(BCS) superfluid to a Bose- Einstein condensation (BEG). 
By contrast, in a Fermi system with p-wave interactions, 
a quantum phase transition is anticipated to take place 
instead of a smooth crossover [l|-|3. 

By sweeping the magnetic field, p-wave Feshbach res- 
onances have recently been observed by several exper- 
imental groups [J-lS|- The binding energies, lifetimes 
and radio-frequency spectroscopy of p-wave Feshbach 
molecules and some low-dimensional aspects of p-wave 
interactions have been investigated in these experiments. 
It is well known that in the dilute, low-energy limit, a 
two-particle s-wave interaction can be well characterized 
using a single parameter, the scattering length. However, 
in p-wave systems a new length scale - the effective range 
of interactions - is required, in addition to the scatter- 
ing volume [l|, [l3 ■ To avoid complicated calculations in 
many-body systems with p-wave interactions, some theo- 
retical descriptions |lll[l3| only use a constant scattering 
volume. This is based on the idea of a pseudopotential, 
which was first introduced by Huang and Yang [13| . Re- 
cently, there have been several developments of the pseu- 
dopotential method for higher-partial wave interactions 
using an energy-dependent scattering length. 
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In particular, numerically exact solutions of two interact- 
ing fermions in a harmonic trap have been studied with 
either s-wave or p-wave couplings. 

Given the interesting and unusual nature of p-wave in- 
teractions, in this paper we aim to investigate the high- 
temperature thermodynamics of a many-body, strongly 
interacting p-wave Fermi gas in a harmonic trap. The 
s-wave high-temperature thermodynamics will be also 
included for completeness. Our investigation is based 
on the virial expansion of the thermodynamic poten- 
tial, which provides a natural bridge between few-body 
and many-body physics |23 - |26| . Restricting ourselves to 
the leading order of the expansion, we calculate the sec- 
ond virial coefficient from two-fermion solutions. Gon- 
sequently, we obtain the high-temperature equation of 
state in the strongly interacting regime, giving the en- 
ergy and entropy as a function of temperature. We ad- 
dress in particular the role of the effective range of p-wave 
interactions in the thermodynamic state equation. 

The calculation of thermodynamics in strongly inter- 
acting quantum many-body systems is always a challenge 
|27H30| . The only conclusive method seems to be that of 
ab — initio quantum Monte Garlo simulation. This, how- 
ever, often suffers from the sign problem for fermions. 
Our virial expansion approach to thermodynamics can 
provide an accurate benchmark for these simulations at 
high temperatures. Already, for s-wave interactions, the 
virial expansion has already been shown to be very useful 
in understanding recent thermodynamic measurements. 
We anticipate that our p-wave virial prediction will also 
give valuable insights and possible calibration procedures 
for future experiments on a Fermi gas near a p-wave res- 
onance. 

The calculation of the second virial coefficient requires 
the whole energy spectrum of the two-fermion bound 
states. To facilitate this calculation, we characterize 



the pseudopotential using two parameters, the scattering 
length and the effective range. This gives a simphfied de- 
scription for an energy-dependent scattering lengtli. We 
validiate our two-parameter pseudopotential treatment 
by comparing the resulting spectrum with the spectrum 
of the full energy-dependent pseudopotential. 

This paper is organized as follows. In the next section, 
we briefly review the derivation of the pseudopotential 
method for all partial-wave interactions and then intro- 
duce the two-parameter pseudopotential. The necessity 
of including a finite range for p-wave interactions is em- 
phasized. In Sec. Ill, we present exact solutions for the 
energy spectrum of two-fermion systems with either s- 
wave or p-wave interactions in a harmonic trap. In Sec. 
IV, we calculate the second virial coefficients. Then, in 
Sec. V we investigate the high-temperature thermody- 
namics of a strongly correlated trapped Fermi gas. Fi- 
nally, Sec. VI is devoted to conclusions and final remarks. 
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tpa (r) = X! X! *-^'™ iJiif^"^) ~ ta.n Sini{kr)] Yi,n{0, ip). 

1=0 m=-l 

(2) 
Here Y;m(0, ip) are spherical harmonics, ji{kr) and ni{kr) 
are spherical Bessel and Neumann functions respectively, 
and Si is the phase shift of Z-th partial wave determined 
by the boundary condition. 

Following Huang and Yang [ij|, we extend Ri (r) — 
ji{kr) — tan Sini{kr) to the vicinity of the origin at r = 0. 
The real interaction potential V^(r) is now replaced by 
the pseudopotential Vps (r). We have 



Vps{r),p{v)^ — {V' + k')^a{r) 



(3) 



The asymptotic functions oi ji{kr) and ni{kr) at r — )■ 
are given by, 



II. TWO-PARAMETER PSEUDOPOTENTIAL 

The first step in understanding the sophisticated 
physics of a quantum many-body system is to model 
the fundamental two-body interactions. At sufficiently 
low temperatures the Heisenberg uncertainty principal 
means that particles must spread over a distance much 
larger than the range of the interaction potential. In this 
case the details of the true interparticle potential become 
irrelevant and it may be modeled by a pseudopotential. 
The basic idea of the pseudopotential is to replace the 
real interaction potential by zero-range potential, which 
acts only at r = and reproduce the same asymptotic be- 
havior of the wave function as would occur with the real 
potential. The first attempt to derive a generalized pseu- 
dopotential was made by Huang and Yang [iJl ■ This has 
been improved by a number of authors. In this section 
we briefiy review the derivation of the pseudopotential 
method 
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A. Two-particle scattering 

We consider a two-particle scattering process and as- 
sume that V{y) is the interaction potential. The motion 
of the center-of-mass can be separated from the relative 
motion, and the Schrodinger equation in the relative co- 
ordinate can be written as 
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respectively. Thus, only ni(kr) is singular at small r. To 
solve Eq.(I3]); we use the method introduced in Idziaszek's 
paper [10] and obtain the pseudopotential for all partial 
waves as 



v,si^{v) = Y.Y. 
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The partial differential operator Yim{d) is obtained 
from the harmonic polynomial r^Yim{f) by replacing the 
Cartesian coordinates Xk with the partial derivatives dx^ 
|3l| . According to Eq. ((4]) , we can easily obtain the form 
of pseudopotential for any specific partial wave. Using 
the expression 



— {y^ + e)^:{v)^V{v)^{v). 
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where fc^ — 2^E/'h? and /i — m/2 is the reduced mass, 
assuming identical particles apart from spin. The wave- 
function has the following asymptotic behavior outside 
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with Ck = {-)''{2l - 2k)l/ {2}k\{l - k)\{l - 2/fc)!), we can 
write the pseudopotential for s- and p-wave as 



and 



27r?i^ tan (So ^, , d 






(7) 



(8) 



respectively. The symbol V (V) acts only to the left 
(right) side of the pseudopotential. 



B. Effective range parameters 

For low energy scattering off a short range potential, 
it is possible to express the variation of the phase shift 
as 
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with only two parameters, the l-th partial wave scattering 
length ai and the effective range r/ [32] . If we define R 
as the finite range of the potential, the effective range ri 
can be expressed as |33| . 
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In the low energy limit we replace the term tan (5;/fc2'+^in 
Eqs. dZl) and ® by using Eq. dH). We then obtain 
the two-parameter pseudopotential for s-wave and p-wave 
interactions as 



Vs = — 
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respectively. 

Neglecting the term rok'^/2 in Eq. IHJ, the two- 
parameter s-wave pseudopotential will result in the fa- 
mous Huang- Yang s-wave zero-range pseudopotential. 
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However, as we know that using the single parameter, 
s-wave scattering length, pseudopotential is not always a 
good approximation in a high-density or tightly trapped 
system, as shown by by Blume and Bolda et al. [20, 34]. 
In such cases, the next term in the expansion needs to be 
included. 

For higher partial- wave interactions, for example p- 
wave scattering, the zero-range pseudopotential fails to 
describe the scattering even in the low-energy limit |10|. 



This can be understood by the zero-range approximation 
constraint condition |l/a^'^^| ^ r;A 
of the effective range equation p^ . 



constraint condition |l/a^'^^| ^ rik'^ or, using the result 



|fca/| < (fci?)2'+i 



(14) 



Writing this out explicitly in the two cases of s-wave and 
p-wave interactions, we see that this implies: 



\kao\ < l/(kR) 
\kai\ < {kRf'^ 



wave case 



p — wave case 



(15) 



This means that zero-range approximation is useful for 
short range interactions with kR — > in the s-wave case, 
since the dimensionless interaction range (kR) can al- 
ways be made arbitrarily small in order to satisfy the 
above inequality, even if \kai\ ^ 1. In practise, this is 
achieved at low density and ultralow temperatures, ie, 
by reducing k. However, this limit cannot be used in the 
strongly interacting p-wave regime with |fca;| ^ 1. For 
these higher partial- wave terms, the right-hand side will 
decrease as kR — > 0, meaning that a finite range cor- 
rection is required in order to reach strongly interacting 
regime. 

It should be noted that for any given values of range 
and scattering length, one can always reach the zero- 
range regime at sufficient dilution, as expected. In the 
s-wave case one has to satisfy k^ <^ {R\ao\)~^ . It is 
clearly possible to reach a regime where simultaneously 
kR <C 1 and |fcao| ^ 1, provided k is small and |ao| is 
very large. However, in the p-wave case the inequality 
becomes 

k^ <t:R/\ai\^ . 

The problem is that if i? ^ |ai|, which is the case near 
a Feshbach resonance, then satisfying this p-wave in- 
equality requires low enough densities such that |A:ai| <C 
|fci?| ' . At such low densities one is no longer in the 
strongly interacting regime, since this inequality is only 
satisfied if \kai\ ^ 1. Finally, it should be pointed out 
that the pseudopotential is not a Hermitian operator. 
However, away from the unphysical region of the origin, 
we have a well-defined scattering problem. 



III. TWO FERMIONS IN A HARMONIC TRAP 

The two-body problem with s-wave and p-wave 
interactions[35] in a three-dimensional harmonic trap was 
solved using a zero-range pseudopotential by Busch and 
Idziaszek et al [3g, [33] . Here we use the two-parameter 
pseudopotential given above. In an axially symmetric 
harmonic trap, the motion of the center of mass can 
be separated from the relative motion, and the relative 
Hamiltonian is given by 



2fi 



'^ 2/22 



') + Vps 



^'. 



rel 



Erel.lpreh 

(16) 



where p = p^ — P21 2: = Zi — Z2 and r = ri — r2 are the 

relative coordinates, /z is the reduced mass, uj is the trans- 
verse frequency of the trap, ry — ijjp/oj, and V^s(r) is the 
two-parameter pseudopotential. To solve Eq. ([T5)) . we 
expand the relative wavefunction ^^ei (r) into the com- 
plete set of the eigenfunctions of three-dimensional har- 
monic oscillator, 
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wherej3„i„2n3(i") = 4>nAilx/d)4'n2{r]y/d)4>n3{z/d), d = 
y^h/ fiu! is the oscillator length in the transverse direction 
and 0n(^) can be expressed by a Hermite polynomial as 



Considering an isotropic harmonic trap, (/3„j^„2„3(r) can 
be replaced by the spherical harmonic function Yi^ (9, (p) 
and the function i?„; (r), written as 



y:>nlm (r) = Rnl (r) Ylm {&, ^) ■ 



(18) 



Using the generalized Laguerre polynomials 
Ln , we can express the function Rni as 

Rnl ^N^i{r/dfexp{-r^/2<f)L^:,+^'^'>{r^/d^). Here 
Nni is the normalization coefhcient. 



Nnoexp (— r^/2d^) L„ (r^/d^) and the relationship be- 

tween the Laguerre polynomials Ln {x) and the conflu- 
ent hypergeometric function C/, 
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Eq. (PO)) can be simplified to 
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where F (x) is the Gamma function. By examining 
the short-range behavior of the confluent hypergeometric 
function, 



T{-v)U[-v,-y]=-^ 
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and taking the approximation fc^ = ifiErei/^ : we obtain 
a secular equation for the relative energy levels: 
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Here the second term on the right hand side of Eq. ([M)) 
shows the effects of the finite-range interaction potential. 
Taking a zero range limit, with ro = 0, Eq. (IM|) gives 
Busch's earlier result [3^. The relative wave functions 
can be found to be: 



A. s-wave interaction in a 3D isotropic trap 



4>rei (r) (X exp ( -^\ r i-v) U l-v, I ^ ) . (25) 



As a preliminary calculation, we consider two fermions 
with unlike spins interacting via s-wave interactions in a 
3D isotropic harmonic trap, with 77 = Wz/w = 1. The 
relative wavefunction can be written in the form of Eq. 
(fT5)) . For s-wave interactions, we only need to keep / — 
in Eq. (fT5|) as / 7^ modes are not affected by the 
interactions. The relative wavefunction of interest can 
then be written as 



%l)rel (r) = ^C„(^„(r), 



(19) 



where (^„ (r) — Rno (r) Yqq {9, ip). We insert this relative 
wavefunction into Eq. (J16p . Following the derivation by 
Busch et al. [3y| the energy level Erei should satisfy the 
following expression. 
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(2n -I- 3/2) hu! and we have de- 
{2i/ + 3/2) hij. By using Rno = 



In Fig. 1, we give the energy spectrum of two fermions 
with s-wave interactions as a function of the dimension- 
less interaction strength d/ao- Here we consider two 
cases: 

(i) zero-range limit (red dotted line) and 

(ii) ro/d = 2R/d = 0.1 (blue solid line). 

It is easy to see that, except for the lowest bound state, 
the two lines differ slightly. The most significant differ- 
ence appears around the unitarity limit with oq — >■ ±cx). 
Thus, the zero-range approximation appear to be sound 
as far as the lowest bound state is concerned. There- 
fore, the scattering length ag can be used to characterize 
the low- energy s-wave scattering, provided that the finite 
range of the interaction potential is smaller than the char- 
acteristic length of the harmonic trap, d. 

In the limiting case of zero scattering length, we find 
that the asymptotic behavior of the energy spectrum for 
the n-th level can be described by. 
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where n = 0,l,2,...isa non-negative integer. 
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Figure 1: (Color online) The energy spectrum of two ferniions 
with s-wave interactions in a 3D isotropic harmonic trap, as 
a function of the dimensionless interaction parameter d/ao at 
two different effective ranges of interactions: ro/d — (red 
dotted lines) and ro/d = 0.1 (blue solid lines). 



Near the uiiitarity limit, the asymptotic behavior of 
the energy spectrum for the n-th level can be written as 
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where 
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B. p-wave interactions in a 3D isotropic trap 

Let us now consider two fermions with p-wave interac- 
tions in a 3D isotropic harmonic trap. We keep only the 
I = 1 term in Eq. (fT5)) . since other I ^ 1 channels are 
not affected by p-wave interactions. We may then write 
the relative wavefunction, 

^rel (r) = Y^ CnlmRnl (r) Yl^ {9, ^) . (28) 



Substituting Eq. (^5)) into Eq. p^ . we obtain 

{ V^:i^ (r) (aydr^) [r3 V<p„i™ (r)] } ,^o 
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where Eni — {2n + 1 + 3/2) huj is the non-interacting 
relative energy and we have written the interacting 
energy eigenvalue as Erei — {2v + 5/2) tku. Using 
T.myL{(^^V)Yim{e,^) = {21 + 1)1 At: and i?„i = 



Nni [r/d) exp {-r^/2d^) U^ {r^jd^), we obtain 
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Using the identity between the Laguerre polynomials 
Ln (r^) and the confluent hypergeometric function U, 

^ LTir'ld^) 1 f 5 r^\ 

^ [2n + 5/2) huj - Erel 2huj ^ "' \ ""' 2' d^ J ' 

(31) 
and the asymptotic behavior of the confluent hyperge- 
ometeric function at r — > 0, 
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the eigenvalue equation Eq. (|5(I)) leads to 
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r(-i^ + i/2) I f dy 1 , , , / 3\ , , 
r(-.-i) — 8-U)+8-(-^-^H''^'^2)-(^^) 

Here we take the same approximation fc^ = 2fiErei/fi^ 
as in the s-wave case. The corresponding eigenfunctions 
have the form. 
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In Fig. 2, we give the energy spectrum of two fermions 
with p-wave interactions as a function of the dimension- 
less inverse scattering volume d^ /a\. The red dotted line 
and blue solid line correspond to a zero effective range 
rid = and a finite effective range rid — —40, respec- 
tively. The value of rid — —40 is calculated from the 
finite range of interaction potentials, R = 0.05d, using 
R = — 2/ri. Unlike the s-wave case, near the unitarity 
limit, a small finite-range parameter in the interaction 
potential can induce a significant change of the energy 
spectrum, implying that the use of zero-range pseudopo- 
tential is not justified for p-wave interactions in the uni- 
tarity limit, as expected. 

This can be understood using Eq. p4p . For p-wave in- 
teractions, the zero-range approximation is valid only if 
|fcai| ^ {kR) ' . For a small value of R and a large scat- 
tering length fli, the wavelength or the collision energy 
therefore should be extremely small in order to satisfy 
this constraint. As explained in Section II, this limit can 
only be achieved if |fcai| <C 1 , which is no longer in the 



C. p-wave interactions and 3D anisotropic traps 
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Figure 2: (Color online) The energy spectrum of two ferniions 
with p-wave interactions in a 3D isotropic harmonic trap, as a 
function of the dimensionless inverse scattering volume d^/af 
at two different effective ranges: rid = (red dotted lines) 
and rid — —40 (blue solid lines). 



Let us now turn to the case of two fermions with p- 
wave interactions in an anisotropic 3D harmonic trap 
/iw^ (rj^p^ + z^) /2. We have to take the general form 
Eq. p7p for the relative wavefunction iprei (r). By sub- 
stituting Eq. (fT7| into Eq. (|16p . we follow the derivation 
for s-wave interactions given by Busch et al. [36|. This 
leads to the following expression for the eigenvalues Erei , 
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(38) 
and the non- 
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interacting energy spectrum is given by En-^n2n3 — 
(ni + 1/2) huj^ + {n2 + 1/2) ^^^ + (^3 + 1/2) hw:,. To pro- 
ceed, we re-write the summation in the above equation 
by using the identity. 



strongly-interacting regime. In contrast, for s-wave in- 
teractions, we only require that |A;ao| <C (fci?) , which 
can generally be satisfied at low dilutions near a Fes- 
hbach resonance. We have checked our two-parameter 
pseudopotential method, by comparing its energy spec- 
trum with that predicted by a more complex energy- 
dependent pseudopotential approach [lOJ. For typical 
parameters corresponding to the experimental condition 
for '^''K fermions, we find an excellent agreement between 
these two pseudopotentials. 

In the non-interacting limit of zero scattering length, 
we find that the following asymptotic form for the n-th 
energy level. 
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In the zero range limit, Eq. (|35p implies that the in- 
teracting energy spectrum differs mostly from the non- 
interacting spectrum at the resonance position oi — > cx). 
However, the existence of a finite range of interactions 
will shift the position to the BCS side, which is deter- 
mined by the condition, (P /{2a\) — rid (n + 5/4) — 0. 

Near the unitarity limit, the lowest energy level is ap- 
proximately zero. The asymptotic value of energy spec- 
trum for the n-th level is given by. 
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ioT E < Eo = h {ujx + ^y + t^z) /2, and use the relation 
satisfied by the Hermite polynomial, 
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As a result, the left-hand side of Eq. ([55|) can be ex- 
pressed by using the function {I — x,y, z ) 
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At r — > 0, the asymptotic behavior of 

Hi {e = {E- Eo) /fuv, r) is given by, 
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and rjx = rjy = rj and rj^ — 1. After a straightforward but 
lengthy calculation, we find that 
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Figure 3: (Color online) The energy spectrum of two ferniions 
with p-wave interactions in a 3D anisotropic harmonic trap, 
as a function of the dimensionless inverse scattering volume 
d /a-i aX ri = Up/uj — 5 and the effective range rid — —40. 



Here we take the same approximation fc^ = 2jj,Erei/fi^ as 
before and Fi {Erei,0) = Fi {EreuY — i)). In an axially 
symmetric trap, we may use the good quantum number 
of projected angular momentum m to label the energy 
levels. The secular equations for ra — ±.1 and for m = 
are, respectively |1C| |. 
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(41) 



(42) 



To calculate the eigenvalue, for Erei < Eq we use the 
integral equation Eq. ([551) . while for Erei > Eq we use 
the recurrence relations [lOj, 



F, (e + 277) - 2F, (e) + F, {e 
^Vl r(-r;/2-e/2) 
2 r(l/2-r;/2-e/2) 



277) 



and 



F,(e + 277,0)-F,(e,0) = -V^ 



r (1/2 -7, -6/2) 



(43) 



(44) 



r (-77-6/2) 

where e = (E — Eq) /hu. 

In Figs. 3 and 4, we show the energy spectrum of two 
interacting fermions with p-wave interactions, as a func- 
tion of the dimensionless inverse scattering volume d^ /a\ 
at an effective range rid = —40, in a cigar-shaped trap 
(77 = 5) and in a pancake-shaped trap (7/ = 0.2), respec- 
tively. These plots are similar to the energy spectrum of 
two fermions in an isotropic trap. However, due to the 
trap anisotropy, the degeneracy for different projected 
angular momentum m is removed. 



IV. VIRIAL EXPANSION OF STRONGLY 
CORRELATED FERMIONS 

The knowledge of few-particle exact solutions provides 
a useful input for investigating the high-temperature be- 
havior and many-body physics of a strongly correlated 



Figure 4: (Color online) The energy spectrum of two fermions 
with p-wave interactions, as in Fig. 3, but with rj = Up/uj = 
0.2. 



quantum gas. This is provided by the quantum virial 
expansion technique [22|-[24l |38J . The essential idea of 
the quantum virial expansion is that at high tempera- 
tures the chemical potential ji is negative and the fugac- 
ity z = exp{fi/kBT) = exp(/3/i) ^ 1 is a well-defined 
small parameter. We can therefore expand the thermo- 
dynamic potential fi of a quantum system in powers of 
the fugacity, regardless of the strength of the interactions. 
In general we may write, 

n = -keTQi {z + 62^2 + • • • + 6„z" + • • • ) (45) 

where Qn is the partition function of a cluster that con- 
tain n particles 



Qn = Tr„ [exp {-U/kBT)] 



(46) 



The trace Tvn takes into account all the n-particle states 
with a proper symmetry. The 7i-th virial coefficient 6„ 
has the form: 



&3 - 



O2 - Q\/2 



Q3-QiQ2 + Ql/^ 
Qi 



(47) 
(48) 



All the other thermodynamic properties can then be de- 
rived from 51 via the standard thermodynamic relations. 
It is convenient to focus on the effect of interactions 
on the virial coefficients. To this end, we consider the 

differences AQ„ = Qn ~ Qn and A6„ = 6„ — &„ , where 
the superscript "1" denotes an ideal, non-interacting sys- 
tem having the same fugacity. Accordingly, we rewrite 
the thermodynamic potential in the form. 



^ = fl^i) - ksTQx [Ab2z'^ + 



A6„z" 



(49) 



where fl^^' is the non-interacting thermodynamic poten- 
tial and 



A62 = AQ2/Q1 
A63 = AQ3/Q1 



Ag2 



(50) 
(51) 



We now describe how to calculate the non-interacting 
thermodynamic potential il^^' and the virial coefficients 

Abr,. 



length and zero effective range of interactions, the spec- 
trum is known exactly: E^ei — (2n -I- 1/2) hu, giving rise 
to, 



A. Ideal thermodynamic potential 

The non-interacting virial coefficients can be deter- 
mined straightforwardly from the non-interacting ther- 
modynamic potential. For a harmonically trapped Fermi 
gas, the non-interacting thermodynamic potential in the 
semiclassical limit {i.e., neglecting the discreteness of the 
spectrum) is, 



C!(i) 



(M 2 7 ^ ^ ' 





(52) 



where Qi = 2 (ksTf / {hwf and hw = rf/^huj. Taylor- 
expanding the non-interacting thermodynamic potential 
in powers of z gives rise to 



h^} 



(-1) 



n+l 



(53) 



B. Second virial coefficient in a harmonic trap 

To calculate the second virial coefficient of a trapped 
interacting Fermi gas, we use Cj — hw /ksT ^ 1 to char- 
acterize the intrinsic length scale relative to the trap. To 
obtain A62, we consider separately AQ2 and Qi. The 
single-particle partition function Qi can be determined 
by the single-particle spectrum of a 3D harmonic oscilla- 
tor or alternatively can be derived from the well-known 
ideal thermodynamic potential. We find 



Qi'^2{kBTf /{hwf 



(54) 



The prefactor of two accounts for the two possible spin 
states of a single fermion. In the calculation of A(32, it is 
easy to see that the summation over the center-of-mass 
energy level gives exactly Qi/2. Thus, the second virial 
coefRcient is determined entirely by the relative energy 
spectrum E^ei and we find that. 



A&2 



\ E hp(- 



■Erei) - eM-E'^^i) 



(55) 



where E^J^ is the relative energy for non-interacting two- 
fermions and the summation is over the whole energy 
spectrum. In the following, we focus on the most inter- 
esting case of the unitary limit. 



A62 = - 



1 cxp (— ci)/2) 



2 [1 + exp (-w)] 



1 
32" 



cD2 



(56) 



The leading term 1/4 on the right-hand side of the 
above equation is universal and temperature indepen- 
dent. This is anticipated following the universality ar- 
gument first suggested by Ho [23|- The term w^ in Eqs. 
([5S| is non-universal and is related to the intrinsic length 
scale of the harmonic trap. However, it is negligibly small 
for a Fermi cloud with a large number of atoms. 

The existence of a finite range interaction term causes 
another non-universal correction. Near resonance, the 
relative energy levels are given by Eq. (I^S|) . 
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where 
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As q ^ 1 for a small effective-range of interactions near 
unitarity, we have 
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(62) 



Thus, the second virial coefficient near resonance can be 
analytically written as 



Ab. 
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kpro' 



32 ^/2^kFao 



'2t: 



(63) 



s-wave interaction 



Let us first consider the simplest case of s-wave inter- 
actions. At resonance with an infinitely large scattering 



where T = T/Tp is the reduced temperature, kp — 

y/2mu}/h{3N)^^^ and Tp ^ {3N)^^^ fuv/ks are Fermi 
momentum and Fermi temperature, respectively. The 
subscript "s" stands for the s-wave interaction. 



2. 



-wave interaction 



For 3D isotropic harmonic traps, the relative energy 
spectrum for near resonant p-wave interactions is de- 
scribed by Eq. 



Erel — 



2n+- + 2q]huj, 



where 



-r(n + 5/2) 



7r(n + 5/4)r(n + l) 



4 
rid 



2d 



ilrf (n + 5/4) 



Note that the lowest energy level is approximately zero. 
In the limits of (rid) — ^ and oi — > oo, the relative en- 
ergy level can be written as Erei = {2n + 5/2) hu, which 
has exactly the same form as the non-interacting spec- 
trum. However, there is a constant shift due to strong 
attractions. As a result, the second interacting energy 
level lines up with the first non-interacting energy level 
and so on. Their contributions to the second virial coef- 
ficient cancel exactly with each other. In the end, only 
the lowest energy level contributes to the coefficient. We 
thus obtain immediately. 



Ab2,p ^^{'21- 



r>^l. 



(64) 



which is universal and temperature independent. The 
factor 21 + 1 accounts for the spectrum degeneracy. 

Near resonance with a small interaction range, we 
have g <C 1. Due to the small size of q, to a good 
approximation, we can use exp [—uj (2n + 5/2 + 2q)] ~ 
(1 -I- 2qu) exp [-w (2n + 5/2)], so that: 
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Then, the second virial coefficient near resonance is given 
by, 
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(65) 

For the 3D anisotropic case, it is difficult to derive the 
asymptotic equation for the second virial coefficient near 
resonance analytically. We have to solve the energy lev- 
els for 7TT, = and m — ±.1 respectively using Eqs. (J4T|) 
and P^ . and calculate the second virial coefficient nu- 
merically. In Fig. 5 we calculate A62.P as a function 
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Figure 5: (Color online) The second virial coefficient A62fp as 
a function of 1/ {kpai) with finite range potential R = 0.05d 
in cigar-shaped trap ry = 5 and pancake-shaped trap r] — 0.2. 



of the dimensionless interaction parameter 1/ (kpai) at 
three different temperatures and a fixed small finite range 
of potential, R = 0.05d. Here we consider a gas with 
N = W^ atoms and use the Fermi temperature as the 
unit for temperature. All the curves with different tem- 
peratures appear to cross at 1/ (kpai) — > 0. This is 
the manifestation of universal behavior anticipated if the 
finite-range corrections are small, meaning that there is 
no intrinsic length scale. 

In Fig. 6 we show the second virial coefficient A62,p 
at unitarity limit oi — )> ±00 as a function of dimen- 
sionless finite range of potential R/d. In the zero-range 
limit A62,p approaches the universal value of 3/2 as 



w = (3iV) 



-1/3 



0.007 is fairly small. 



V. HIGH-T THERMODYNAMICS OF 
STRONGLY CORRELATED FERMIONS 

We are now in position to study the equation of state in 
the high-temperature regime. Using the thermodynamic 
relation, we have the number equation, 



iV = 7V(^' +2 



{ksTf 



and the total energy. 



(2A62Z^ + 3A6, 



3Z 



E ■ 



{hwf Tp' ' 



(66) 



(67) 



where A6„ = dAb„/dT and the non-interacting number. 



iV(i) = - 



2 jkBTy 1 
{hwf 2 
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(68) 
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Figure 6: (Color online) The second virial coefficient A62,p as 
a function of the dimensionless finite range Rjd in the uni- 
tarity limit in a cigar-shaped trap (r; = 5) and in a pancake- 
shaped trap {r\ = 0.2). 



Figure 7: (Color online) Temperature dependence of the en- 
ergy of a strongly correlated Fermi gas with s-wave and p-wave 
interactions. For comparison, we plot the ideal gas result us- 
ing solid lines. 



The entropy of the system can be calculated using, 
S={E-^- fj.N) /T, 



(69) 



where the chemical potential fi = fesTlnz (not to be 
confused with the reduced mass). Eqs. (|52|) . ([GG)) . (j67)) . 
together with ([68l) . form a closed set of expressions for 
thermodynamics, which can be solved self-consistently. 

Up to the second order virial expansion, the energy 
and entropy of a Fermi gas with resonant s-wave and 
p-wave interactions are shown in Figs. 7 and 8, re- 
spectively. Here, we consider a zero-range potential so 
the thermodynamics is universal. The Fermi energy is 
Ep = kgTp ~ {"iN) ' tujj. For comparison, we show 
the ideal gas result using the solid lines. It is clear that 
the equation of state are strongly affected by interactions, 
even in the high temperature regime. This interaction ef- 
fect is particularly apparent for p-wave interactions due 
to its large second virial coefficient. 

In Figs. 9 and 10, we present the non-universal effect 
on thermodynamics caused by the presence of a finite- 
range of interaction potentials. Eqs. (|63l) and (|55t are 
used in the calculations of the virial coefficient with the 
inclusion of the influence of the flnite-range of interac- 
tions. 



VI. CONCLUSIONS AND REMARKS 

In conclusion, we have investigated the high- 
temperature thermodynamics of a strongly correlated 
Fermi gas with s-wave and p-wave interactions, by ap- 
plying a second-order quantum virial expansion method. 
The second virial coefficient has been calculated, based 
on the numerically exact solution for two fermions in a 
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Figure 8: (Color online) Temperature dependence of the en- 
tropy of a strongly correlated Fermi gas with s-wave and p- 
wave interactions. The ideal gas results are shown by the solid 
lines. 



harmonic trap. In this study, we have particularly fo- 
cused on the effects arising from a finite-range interaction 
potential, which is crucial for p-wave interactions near the 
unitarity limit. Non- universal corrections due to the fi- 
nite range corrections to the thermodynamics have been 
addressed. 

We expect that these thermodynamic results should be 
useful as a complementary approach to quantum Monte 
Carlo simulations in understanding experimental ther- 
modynamic measurements for a p-wave Fermi gas near a 
Feshbach resonance. 
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